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loaded into a one-dimensional optical lattice are investigated by means of a conformal field 
theory approach. We show that for attractive interactions two different superfluid phases 
emerge fori^>3/2: ABCS pairing phase, and a molecular superfluid phase which is formed 
from bound-states made of 2^^^ fermions. In the low-energy approach, the competition be- 
tween these instabilities and charge-density waves is described in terms of Z^v parafermionic 
degrees of freedom. The quantum phase transition for F = 3/2, 5/2 is universal and shown 
to belong to the Ising and three-state Potts universality classes respectively. In contrast, for 
F > 7/2, the transition is non-universal. For a filling of one atom per site, a Mott transition 
occurs and the nature of the possible Mott-insulating phases are determined. 
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I. INTRODUCTION 

Ultracold atomic physics has attracted a lot of interest in recent years with the opportunity 
to study strongly correlated effects, such as high-temperature superconductivity, in a new context 
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In particular, loading cold atomic gases into an optical lattice allows for the realization_of 
bosonic and fermionic lattice models and the experimental study of exotic quantum phases 
A prominent example is the toolbox to engineer the three-dimensional Bose-Hubbard model 
and the observation of its Mott insulator - superfluid quantum phase transition with cold bosonic 
atoms in an optical lattice ^. 

Ultracold atomic systems also offer an opportunity to investigate the effect of spin degeneracy 
since the atomic total angular momentum F, which includes both electron and nuclear spins, can 
be larger than 1/2 resulting in 2F-\-l hyperfine states. In magnetic traps, these 2F + 1 components 
are split, while in optical traps these (hyperfine) spin degrees of freedom are degenerate and novel 
interesting phases might emerge. For instance, Bose-Einstein condensates of bosonic atoms with 
a nonzero total spin F are expected to display rich interesting structures in spin space like in 
superfluid ^He. In this respect, various exotic superfluid condensates (including nematic ones), 
Mott insulating phases, and_non-^rivial_ vortex structures have been predicted recently in spinor 
bosonic atoms with F > 
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These theoretical predictions might be checked 



in the context of Bose-Einstein condensates of sodium, rubidium atoms [12] and in spin-3 atom 

ike ^Li, '^'^K or ^'^^Yb, is also expected 



of Cr [13]. The spin-degeneracy in fermionic atoms^ 
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In particular, a molecular 



to give rise to some interesting superfluid phases 
superfluid (MS) phase might be stabilized in multicomponent attractive Fermi gas where more 
than two fermions form a bound state. Such a non-trivial superfluid behavior has already been 
found in different contexts. In nuclear physics, a four-particle condensate like the a particle is 
favored over deuteron condensation at low density J[19|] and it may have implications for light 



nuclei and asymmetric matter in nuclear stars 
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lis quartet condensation can occur in 



the field of semiconductors with the formation of biexciton 



22l |. A quartetting phase, which stems 



:'rom the pairing of Cooper pairs, also appears in a model of one-dimensional Josephson junctions 



2^. A po.*le expenmental „b„™ti„n of ,u„tets .ight be found in snpo.condnctn.g ^uantnn, 

interference devices with (100)/(110) interfaces of two d-wave superconductors [2j|. In particular, 
the he/ 4:6 periodicity of the critical current with applied magnetic flux has been interpreted as 



the formation of quartets with total c harg e 4e 



predicted in four-leg Hubbard ladders 
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251]. Finally, a bipairing mechanism has also been 



Recently, the emergence of quartet and triplet (three-bo 
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)ound states) has been proposed to 
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32t | . Much of these studies 



occur in the context of ultracold fermionic atoms 
have been restricted to the special case of an SU(n) symmetry between the hyperfine states of a 
n-component Fermi gas. In this paper, we investigate the generic physical features of half-integer 
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spins F = N — 1/2 fermionic cold atoms with s-wave scattering interactions loaded into a one- 
dimensional optical lattice. The low-energy physical properties of 2F+1 = 2A^-component fermions 
with contact interactions are known to be described by a Hubbard-like Hamiltonian Qj: 

J 

^ = E + H.c] -f^Y.'^i + Y.Uj PjM,iPjM,i, (1) 

i,a i i,J M=—J 

where ^ (a = 1,...,2A^) is the fermion creation operator corresponding to the 2F + 1 = 2N 
atomic states and Ui = Yla is the density operator on site i. The pairing operators in 

Eq. ([T|) are defined through the Clebsch-Gordan coefficients for forming a total spin J from two 
spin-F fermions: Pj^i = Ylaf3(-^-^\-^^-^'^^P)^^ai'^\i- '^^^ interactions are SU(2) spin-conserving 
and depend on Uj parameters corresponding to the total spin J which takes only even integers 
value due to Pauli's principle: J = 0,2, ...,2N — 2. Even in this simple scheme the interaction 
pattern is still involved since there are coupling constants in the general spin-F case. In this 
paper, we shall consider a two coupling-constant version of model ([T]) with U2 = ■■■ = U2N-2 7^ Uq 
which incorporates the relevant physics of higher-spin degeneracy with respect to the formation of 
an exotic MS phase. To this end, it is enlightening to express this model in terms of the density 
Ui and the BCS singlet-pairing operator for spin-F which is defined by: 



a, 13 ^ a 

where the matrix J is the natural generalization of the familiar antisymmetric tensor e = ia2 to 



spin F > 1/2. Such a singlet-pairing operator has 
dimensional frustrated spin- 1/2 quantum magnets 



jeen extensively studied in the context of two- 
Using the relation Y.j,m ^jm i^JM,i = nf. 



model ([I]) with U2 = ■■■ = U2N-2 / Uq reads then as follows: 

i,a ill 

with U = 2U2, V = Uo-U2, and fi = fi- U/2. 

When V = 0, i.e. Uq = U2, model ^ corresponds to the Hubbard model with an SU(2A^) spin 
symmetry: Ca,i — > X]/3 ^apcp^i, U being a SU(2A^) matrix. The Hamiltonian (I2D for F / still 
displays an extended symmetry. Indeed, since the BCS singlet-pairing operator ([2|) is invariant 
under the Sp(2A^) group which consists of unitary matrices U that satisfy U*J'W = J , model 
^ acquires an Sp(2A^) symmetry [3j]. In the spin F = 1/2 case, model ([HD reduces to the SU(2) 
Hubbard chain since SU(2) ~ Sp(2). In the F = 3/2 case, i.e. N = 2, there is no fine-tuning 
and models ([1]) and ([3]) have thus an exact Sp(4)~ S0(5) symmetry [35|. For F > 3/2, the 
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Sp(2A^) symmetry is only present for the two coupling-constant model ([3]). The Sp(2A^) symmetry 
simplifies the problem but may appear rather artificial at this point. However, one reason to 
consider this model stems from the structure of the interaction involving the density rij and the 
BCS singlet-pairing term ([2]) which are independent for F > 3/2 ^. The competition between these 
two operators gives rise to interesting physics. Indeed, at [/ = and large V < 0, a, quasi-long- 
range BCS phase, which is a singlet under Sp(2A^), should emerge. In contrast, for V = and 
U < 0, we expect the stabilization of a charge-density wave (CDW) or a MS instability made of 
2N fermions Mi = c\-. . . c^2N i which is an SU(2A^) singlet. We thus expect the existence of several 
distinct phases in the Sp(2A^) model ([3]). As it will be revealed in the following, the competition 
between these instabilities relies on the existence of a Zjv symmetry which is also a symmetry of 
the original SU(2) model ([T|). This Z^v symmetry can be defined by considering the coset between 
the center 'L2N of the SU(2A^) group ^ and the center Z2 of the Sp(2A^) or SU(2): Z^r = 1i2n /'^2 
with 

'^2N : - e^""/^cl,, , n = 0, 2iV - 1 (4) 
Zjv:<i - e^"-/^4^,,n = 0,....,iV-l, (5) 

the Z2 symmetry being Ca,i —Ca,i- The Zjv symmetry ([5]) provides an important physical 
ingredient which is not present in the F = 1/2 case and stems from the higher-spin degeneracy. The 
stabilization of a quasi-long-range BCS phase for F > 1/2 requires the spontaneous breaking of this 
Zat symmetry since the singlet pairing Pqq^ ^ is not invariant under this symmetry. In contrast, 
if Ztv is not broken, a MS phase, which is a singlet under the Z^r symmetry, might emerge. In the 
following, the delicate competition between these superfluid instabilities will be investigated with 
a special emphasis on this Z^r symmetry. In particular, in the low-energy approach, this discrete 
symmetry is described by the Zjv parafermionic conformal field theory (CFT) which captures the 



universal properties of two-dimensional Z^r generalized Ising models j36l ]. This approach will enable 
us to determine the main features of the zero temperature phase diagram of model ([3|). A brief 
summary of our main results have already been published as a letter [29i] . 

The rest of the paper is organized as follows. The low-energy effective theory corresponding 
to model ^ is developed in section II. The nature of the different phases is then deduced from a 
renormalization group (RG) analysis. In section III, we present the phase diagram of the model 



^ For F = 1/2, one has indeed PJq i^bo.i = nf — Ui. 
We recall that the center of a group G is defined from the elements which commute with all elements of G. In 
particular, the Sp(2Ai') or SU(2) groups share the same center which is the Z2 group. 
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([3]) at zero temperature for incommensurate filling and we discuss the main physical properties of 
the phases as well as the nature of the quantum phase transitions. The Mott-insulating phases for 
the case of one atom per site are also investigated. Finally, we conclude in section IV and some 
technical details are presented in three appendices. 



II. LOW-ENERGY APPROACH 



In this section, we present the low-energy description of the spin-F Hubbard model with a 
Sp(2A^) symmetry ([3]) which will enable us to determine its phase diagram at zero temperature in 
the next section. 



A. Continuum limit 



Let us first discuss the continuum limit of the lattice model ([3]). Its low-energy effective field 
theory can be derived from the continuum description of the lattice fermionic operators C n^j in 



terms of right and left-moving Dirac fermions (for a review, see for instance the books 371. ISSl. 13911): 



Ca,i/V^ ^ Rae""^" + Lae-"'^^ (6) 

with X = iao (ag being the lattice spacing) and kp is the Fermi momentum. The left (right)-moving 
fermions are holomorphic (antiholomorphic) fields of the complex coordinate z = vpT + ix (r being 
the imaginary time and vp is the Fermi velocity): Ra{z),La{z). These fields obey the following 
operator product expansion (OPE): 



2lT [z — u>) 



In this continuum limit, the non- interacting part of the Hamiltonian ([3]) corresponds to the Hamil- 
tonian density of 2N free relativistic massless fermions: 

Tio = -ivF (^: Rld^Ra ■ - ■ L^dxla :) , (8) 

where vp = 2taQs'm(kpa()) and the normal ordering :: with respect to the Fermi sea is assumed 
as well as a summation over repeated indices. The continuous symmetry of the non-interacting 
part of the Hamiltonian ^ is enlarged, in the continuum limit, to \J{2N)\l ® \]{2N)\ji under 
independent unitary transformations on the Dirac fermions. 
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The next step of the approach is to use the non-Abehan bosonization [40|, l4ll . |42| | to investigate 
the effect of the Sp(2A^) symmetry of the lattice model. To this end, we introduce a SU(2A^)i 
Wess-Zumino-Novikov-Witten (WZNW) primary field g^/s to represent the spin degrees of freedom 
and an additional U(l) charge boson field <I>c- The free-fermion theory ([8]) is then equivalent to 
the CFT with an U(l) x SU(2A^)i symmetry. The right and left-moving fermions can be written 
as a product of spin and charge operators: 



Ra ~ : exp (i\/27r/iV <I>c_r,) : QaR 



Ll ~ : exp (i^2n/N ^^A : 9aL, (9) 



where gu and are respectively the right and left parts of the SU(2A^)i primary field g with 
scaling dimension (2A^ — 1)/2N which transforms in the fundamental representation of SU(2A^). 
In Eq. ([9]), ^cR,L denote the chiral parts of the charge boson field: = ^cR + ^cL- As is well 
known, the free-Hamiltonian ^ can be decomposed into charge and spin pieces: 



{d.^c)' + [d^Qcf] + [: lili ■ + : IM :] , (10) 

where 0c = <^cL — ^cR is the dual charge boson field and I^j^{A = 1, . . . , 4A^^ — 1) are the currents 
which generate the SU(2A^)i conformal symmetry with central charge c = 2N — 1. They express 
simply as bilinears of the Dirac fermions: 

Ir ='■ RaTaP^P ='■ L^a^apl^li (11) 

T"^ being the generators of SU(2A'^) in the fundamental representation normalized such that: 
Tr(T'^T^) = /2 (see Appendix A). These currents satisfy the SU(2A^)i current algebra 

^AB A fABC 
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xAB • fABC 

4 {z) li (o) ~ ^ ' + J;, ^-. /g {^) , (12) 

Svr^ {z - uj) 27r [z - lv) 

where f^^^ are the structure constants of the SU(2A^) group. 

In the presence of interactions, the spin-charge separation (|lUp still holds away from half-filling 
(i.e. N atoms per site). In that case, the low-energy Hamiltonian of model ([3]) separates into two 
commuting charge and spin pieces: 

n = nc + ns, [nc,ns] = o. (is) 

Let us first consider the charge degrees of freedom which are described by the bosonic field $c- 
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1. Charge degrees of freedom 



Using the continuum description ([6]) and the decomposition ([2]), we find that the charge exci- 
tations of model ([3]) are captured by the following Hamiltonian: 



He 



2V + UN{2N -I) ^ ,2 



It can be recast into the form of the Tomonaga-Luttinger free Hamiltonian: 



-He 



Vr. 



(14) 



(15) 



where the charge velocity Vc and the Luttinger parameter Kc read as follows: 

= VF[l + aQ{2V + UN{2N -1))/{Nttvf)]^^'^ 



Kc = [l + aoi2V + UN{2N -l))/{NTrvF)] 



-1/2 



(16) 



The conserved quantities in this U(l) charge sector are the total particle number M and current 
I: 



j dx : RlRa + LILo, : = yj2N/^ j 



dx dr^r 



I = j dx : L^La, - RlRa ■= \/2NJ^ j dx (9^9^. 



(17) 



For generic fillings, no umklapp terms appear and the charge degrees of freedom display metallic 
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391 ]. However, as it will be 



(gapless) properties in the Luttinger liquid universality class 
discussed in section HI, the existence of an umklapp process, for a commensurate filling of one 
atom per site, is responsible for the formation of a charge gap and the emergence of different Mott 
insulating phases. 



2. Spin-degrees of freedom 

All non-trivial physics corresponding to the spin degeneracy is encoded in the spin part TCs of 
the Hamiltonian (jl3p . Its continuum expression can be obtained by decomposing the SU(2A^)i 
currents of Eq. (fTTj) into || and _L parts 1^^^ = {I^^,I\_)r,l with respect to the Sp(2A^) symmetry 
of the lattice model ([3]). The currents IIr{^l)^^ = '^,-,N{2N + 1) generate the Sp(2A^)i CFT 
symmetry with central charge c = N{2N + 1)/{N + 2). They can be simply expressed in terms of 
the chiral Dirac fermions: 



8 



where T'^ are the generators of Sp(2A^) in the fundamental representation and normahzed such 
that: Tt:{T"'T^) = (5*^^/2 (see Appendix A). These currents verify the Sp(2iV)i Kac-Moody algebra 
given by Eqs. (112p with /"''^ the Sp(2A^) structure constants. The remaining SU(2A^)i currents 
L = 1) ■•5 2A^^ — — 1) are written as: 

IlR=:RinpRp:, IlL=:Lin^Lp:. (19) 

With these definitions, we can now derive the continuum description of the spin degrees of 
freedom using Eq. ^ and Eqs. (|A2t lASP of Appendix A. After some cumbersome calculations, 
we find that the low-energy Hamiltonian in the spin sector, 7is, can be expressed in terms of the 
currents only, and displays a marginal current-current interaction: 



~ 2A^ + 1 



ja ja .1 ja ja 



+ [■ IIrAr ■■ + ■■ AlIYl ■■] +9\\I^r%+9±IYrIYl, 

(20) 

with = — 2ao(2y + NU)/N, g± = 2ao(2y — NU)/N and are the spin velocities: Vg\\ = 

vp — UaQ/2-K — VaQ/-K, Vs± = vp — UaQ/2-K + VaQ{N + 1)/A^7r. In the simplest = 1 case, we have 
Sp(2)i ~ SU(2)i and the Hamiltonian (I20p with ^j. = describes the spin-sector of the continuum 
limit of the spin-1/2 Hubbard chain which can be found for instance in the book 3^ written with 
the same notations used in this work. In the N = 2, i.e. the F = 3/2 case, one can express 
model ()20p in a more transparent basis. Indeed, there is a simple free-field representation of the 
unperturbed SU(4)i ~ S0(6)i CFT in terms of six real (Majorana) fermions which has been used 
in the context of the two- leg spin-1/2 ladder with four-spin exchange interactions Introducing 
six real fermions and CRi,i = 1, 5 to describe respectively the Z2, i.e. Ising, and S0(5)i ~ 
Sp(4)i CFTs, the currents of Eq. (I20p can be written locally in terms of these fermions: 

-^±R,L ~ ~-^^R,L^R,L 

im^L - -^^kdkL^ (21) 

where a = 1, . . . , 10 = where 1 < i < j < 5. The Majorana fermions are normalized as the 
Dirac fermions of Eq. ([7]) to reproduce faithfully the Kac-Moody algebra (jl2p . The current-current 
model of Eq. (j20p can then be expressed in terms of these real fermions: 



-Hs = -f [: ■■ - : ClQ^Cl ■] [■ ^R^-^R ■ " ■ ^'9.^1 ^]+\\mif + ^±eRCURa, (22) 

with A|| = —ao{U + V), X± = aQ{V — U), and the spin velocities: v = vp — clo{U + V)/2tt, 
vq = vp — ao{U — 3V)/2Tr. In absence of the spin-velocity anisotropy, model (p2|) turns out to be 
exactly solvable and has also been studied in the context of a SO (5) symmetric two- leg ladder j45l |. 
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B. RG analysis 

In the general N case, model (I20p is not integrable and the main effect of the current-current 
interaction can be elucidated by means of a RG analysis. The one-loop RG equations of model 
([201) are given by (see Appendix B): 

N+1 2 N -1 2 

g± = ^9±9\\, (23) 

where g'^ y = dg±^\\/dt, t being the RG time. In Eq. ([23]) . we have neglected the spin-velocity 
anisotropy Vs = Vs\\ — Vsi_ and have absorbed Vs in a redefinition of the coupling constants: 
Oa 9a /vs- We have also obtained the two-loop RG equations of model (pOj) and the results are 
presented in Appendix B. The one-loop RG equations (j23p can be solved and in particular the RG 
invariant flow K which parametrizes the RG lines reads as follows: 

-{N+l)/N / 2 2 



K = gl'''^'"" [9t-9i)- (24) 

The RG flow emerging from Eqs. (j23p is rich and consists of three different phases (see Fig. 1). 

In region I, where both U and V are positive, all couplings go to zero in the infrared (IR) limit 
and the interaction is marginal irrelevant. The symmetry of the IR fixed point is SU(2A^)i (up to a 



velocity anisotropy) leading to 2A^ — 1 gapless spin excitations. T 



le current-current interaction of 



46l . I47l | . The low-energy properties 



Eq. (j20p leads to logarithmic corrections to physical quantities 
of this phase are very similar to that of the repulsive SU(2A^) Hubbard chain which have been 



studied in Refs. 



49|. In contrast, a spin gap is opened by the interaction in the two remaining 
phases. In phase II, defined by C/ < and V > NU /2, the RG flow in the far IR limit is attracted 
along a special symmetric ray = g± = g* > where the interacting part of the Hamiltonian 
(|20p can be rewritten in a manifest SU(2A^) invariant form: 

Km = 9* + liRlii) = 9*Ipt- (25) 

This IR Hamiltonian thus takes the form of the SU(2A^) Gross-Neveu (GN) model [50,] which is an 
integrable massive field theory [sil. The development of the strong-coupling regime in the SU(2A^) 
GN model leads to the generation of a spin gap. The low-energy properties of the spin sector of 
phase II can be extracted from the integrability of the SU(2A^) GN model ()25p . Its low-energy 
spectrum consists into 2A^ — 1 branches with masses: = m sin(7rr/2A^) that transform in the 



SU(2A^) representation with Young tableau with one column and r boxes (r = 1, . . . , 2A^ — 1) 



5ii. 
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FIG. 1: Phase diagram of model ([20| obtained from the RG approach. Phase I is gapless while phases II 
and III have gapped spin excitations. The dash lines stand for the phase transition lines predicted from the 
RG calculation. 

The corresponding eigenstates are labelled by quantum numbers associated with the conserved 
quantities of the SU(2A^) low-energy symmetry (the Cartan basis): = J dx + -fy'/j), a = 
(1, A^) and Q\ = j dx (/^^ + i = (1, ...,N - 1). Due to the Sp(2A^) symmetry of model 

([3]), the Qjj numbers are conserved whereas the Q\ charges are only good quantum numbers at low 
energy. This is an example of a dynamical symmetry enlargement which corresponds to a situation 
where a Hamiltonian is attracted under a RG flow to a manifold possessing a higher symmetry 
than indicated by the original microscopic theory. This phenomenon occurs in a large variety of 
models with marginal interactions in the scaling limit as, for instance, in the half-filled two-leg 



Hubbard model [53|] and in the SU(4) Hubbard model at half filling [5^, where a SO (8) symmetry 
emerges at low energy. 

In the second spin-gapped phase (HI) of Fig. 1, defined by F < and V < NU/2, the RG flow 
is now attracted along the asymptote: = —g± = g* > 0. In that case, the interacting part of 
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the IR Hamiltonian becomes 

K,int=5*(/^H^|'[L-/lij/iL), (26) 

which can be recast as a SU(2A^) GN model (I25p by means of a duality transformation T> on the 
fermions: 'DR{L)'D-'^ = R{L) with 

Roi = Cfoipl^fji La = (27) 

Using Eqs. p8l [T9|) . we observe that this transformation acts on the currents as: = ^\\r{l) 

and = ~(+)-^1r(l) ^° indeed maps ([26|) onto (p5]) . Besides the opening of a spectral 

gap, we thus find that phase III possesses a hidden symmetry at low energy i.e. a SU{2N) 
symmetry generated by the dual currents {I^jk^^^ , ^^pk^l) ) • "^^^ ^P™ spectrum in phase III can be 
obtained from the duality symmetry D and consists into the 2N — 1 branches rrir which transform 
in the representations of the dual group SU{2N). The dual quantum numbers are now given by: 
Q|| = Q|[ and Q\ = J dx {P^l + Ar) = I dx =T\. We thus observe that the low-lying 

excitations in phase III carry quantized spin currents in the "_L" direction. In this sense, the phase 
III might be viewed as a partially spin-superfluid phase. In summary, the existence of these two 
distinct spin-gapped phases is a non-trivial consequence of higher-spin degeneracy and does not 
occur in the F = 1/2 case. 



C. Conformal embedding 

The crucial point of a non-perturbative analysis is often the identification of a good basis that 
describe low-lying excitations of the phase. Much insight on this problem can be gained from the 
symmetries of the model and the use of the non-Abelian bosonization approach for ID systems. 
Such an approach has 



and in spin chains 
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3een extremely powerful in the past as in quantum impurity problems 



55| 



42]. So far, we have used an U(2iV)i = U(l) x SU(2iV)i CFT approach to 
determine the low-energy properties of model ([3]). However, this description is not adequate to give 
a full understanding of the two spin-gapped phases found in the RG approach. In particular, the 
physical origin of the formation of the spin gaps is not clear at this point. What is the nature of 
the discrete symmetry which is spontaneously broken in phases II and III? A second weak point of 
the previous analysis is the determination of the quantum phase transition between the two spin- 
gapped phases which is rather unclear within the preceding description. At least, this transition 



should occur in the manifold which is invariant under the duality V symmetry (I27p i.e. the self-dual 
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manifold defined by: g±^ = 0. The low-energy Hamiltonian which describes the phase transition is 
thus (neglecting the spin- velocity anisotropy): 

2ttv^ 



n 



SD 



. jci ja , , , ja ja , , , T% ji . \ . 
• -'lIR-'lli? • + • -'IIL-'IIL • + • ^±R^±R ■ + ■ ^±L^±L 



+ 9rVk^ (28) 



2N + 1 

with > so that the marginally relevant current-current interaction opens a mass gap in the 
Sp(2A^) sector. However, one cannot conclude on the occurrence of the first-order phase transition 
since the non-interacting Hamiltonian of Eq. (|28p contains more degrees of freedom than the 
Sp(2A^) one. They remain massless and control the quantum phase transition. The nature of these 
decoupled degrees of freedom is not clear at this point. 

Prompted by all these questions, it is important to fully exploit the existence of the Sp(2A^) 
symmetry of the lattice model ([3]) and to consider the following conformal embedding: U(2A^)i 
^ U(l) X Sp(2iV)i X [SU(2iV)i/Sp(2iV)i]. The coset SU(2iV)i/Sp(2iV)i CFT has central charge 
c = 27V - 1 - iV(2iV + l)/(iV + 2) = 2{N - 1)/{N + 2) which is that of the Zn parafermionic 
CFT [36]. In fact, as shown by Altschuler [sel, it turns out that the SU(2N)i/Sp(2N)i CFT is 
indeed equivalent to the Z^v parafermionic CFT which describes self-dual critical points of two- 
dimensional Zjv Ising models (see Appendix C). The conformal embedding U(2A^)i U(l) x 
Sp(2A^)i X Zjv provides us with a non-perturbative basis to express any physical operator in terms 
of its charge and spin degrees of freedom which are described respectively by the U(l) and Sp(2A^)i 
X Zjv CFTs. Since Sp(2) ~ SU(2), this basis for = 1 accounts for the well-known spin-charge 



separation which is the hallmark of ID spin-1/2 electronic systems [38l. l39l| . In this respect, in the 
general half-odd integer spin case, the Ztv symmetry plays its trick and provides a new important 
ingredient not present for F = 1/2. In the low-energy approach, the spin degrees of freedom 
corresponding to this symmetry are captured by an effective 2D Z^v Ising model. As in the N = 2 
case, these Zjv Ising models exhibit two gapped phases described by order and disorder parameters 
(Tfc and fik (k = 1, ••, — 1) which are dual to each other by means of the Kramers- Wannier (KW) 
duality symmetry. This duality transformation maps the Z^v symmetry, spontaneously broken in 
the low-temperature phase ((ufc) ^ and (fik) = 0), onto a Zjv symmetry which is broken in the 
high-temperature phase where (/i^) ^ and {a^) = 0. At the critical point, the theory is self-dual 
with a "Zn X Z^r symmetry and its universal properties are captured by the Zjv parafermionic 



CFT with (Tk^lJ-k becoming conformal fields with scaling dimension dk = k{N — k)/N{N + 2) [36|]. 
This Zjv CFT is generated by right and left parafermionic currents ^kR,L {^\rl ~ ^N-kR,L, 
k = 1, . . . , N — 1) with scaling dimension = A;(A^ — k)/N which are the generalization of the 
Majorana fermions of the Z2 Ising model. Under the Z^r x Z^v symmetry, ^kL (respectively ^kR) 
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carries a {k, k) (respectively {k, —k)) charge which means: 

^kL,R - e±^"/iv^^^_^ ^^^^^ (29) 

with m = 0, . . . , — 1. As it is discussed in Appendix C, there is a faithful representation of 
the first parafermionic currents ^il,r in terms of the Dirac fermions and the charge bosonic field 

LiJapL^ ~ ^ : exp (V8vr/iV : 

RiJapRl - ^ : exp (-iVSvr/iV : (30) 

We deduce, from this identification, that the duality symmetry (j27p of the RG analysis together 
with the Gaussian duality {^cR,l —>' T^cR,l) give — > and — > — which is nothing 
but the KW duality transformation on the first parafermionic current (see Eq. ()C2p of Appendix 
C). It is thus tempting to interpret the formation of the spin-gapped phases of Fig. 1 as the 
result of the spontaneous breaking of the and Zn discrete symmetries. In this respect, from 
the correspondence (|30]l and Eq. (f29l) . we observe that under the Zat symmetry, the left and 
right-moving fermions transform as: 

L„^e-^™/^L„, Ra, ^ e-''''^/^ Ra, (31) 

while under Ztv we have: 

L„^e-^™/^L„, Ra ^ e'''"'/^ Ra. (32) 

Using the continuum representation of the fermions ([6|), we find that the lattice Z^r symmetry ([5]) 
corresponds to the Ztv symmetry of an effective 2D Zjv Ising model. In contrast, the Z^v symmetry 
has no simple local lattice representation. 

We need now to fully identify phases II and III of Fig. 1 with the high and low-temperature 
phases of the Z^v Ising model. To this end, we use the fact that the SU(2A^)i WZNW primary field 
g of Eq. Q can be expressed in terms of the Sp(2A^)i x Z^v basis using the conformal embedding 
(see Eq. of Appendix C): 

: exp (-VW^^^c) : L^R^ ~ Tr 5 ~ /iiTr (^^^\ (33) 

where (p^^^ is the Sp(2A^)i primary field with scaling dimension A = (2A^-|-l)/2(A^-|-2). We deduce, 
from the identification (|33p . that Tr g has the same behavior as ni under the I^isr >i "^n symmetry 
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(Eq. del]) of Appendix C): 

TV (7 —> TV (7 under Zjv 

Trg^ gi2™/AfTr g under Z^. (34) 

The next step of the approach is to note that the interacting part of the Hamiltonian (j25|) . which 
controls the strong couphng behavior of the RG flow in phase II, can be expressed in terms of 
Tr g: HI jj^^ ~ |Tr g\^ . The ground state of this phase displays long-range order associated with 
the order parameter Tr g: (Tr g) ^ 0. According to Eq. we then deduce that the Z^r 

symmetry is spontaneously broken while the Zjv symmetry remains unbroken in phase II. The Z^r 
Ising model thus belongs to its high-temperature phase and a spectral gap is formed. Using the 
KW duality symmetry or the transformation (j27p on the Dirac fermions, one can conclude that 
phase III corresponds to the low-temperature phase of the Z^v Ising model where the Z^v symmetry 
is spontaneously broken. 

In summary, the existence of the two spin-gapped phases of Fig. 1 is a non-trivial consequence 
of higher-spin degeneracy and does not occur in the F = 1/2 case. The emergence of the spin-gap 
stems from the spontaneous breakdown of the Zjv or Z^v discrete symmetries. As we shall see now, 
these symmetries are central to the striking physical properties displayed by these phases. 

III. PHASE DIAGRAM 

In this section, we discuss the phase diagram at zero temperature of the lattice model ^ for 
incommensurate filling and for a commensurate filling of one atom per site. The nature of the 
quantum phase transitions will also be investigated. Let us start with the incommensurate filling 
case. 

A. Phase diagram for incommensurate filling 

We shall now determine the nature of the dominant electronic instabilities of the different phases 
of Fig. 1 for incommensurate filling. The RG analysis of the preceding section reveals the existence 
of three different phases. 



15 




FIG. 2: Phase diagram of model ^ for incommensurate filling. 
1. Critical phase 

In phase I of Fig. 2, the interaction in the spin sector (Eq. (j20p ) is marginally irrelevant when 
U and V are both positive so that the phase displays an extended quantum critical behavior. 
Up to a spin-velocity anisotropy, the low-energy properties of this phase are very similar to that 
of the repulsive SU(2A^) Hubbard chain with one gapless charge mode and 2N — 1 gapless spin 
excitations Q, [3]. All correlation functions in this phase display a power-law behavior. The 
leading instabilities are those which have the slowest decaying correlations at long distance. For 
phase I, the dominant instabilities are the 2kF CDW p2kp and generalized 2kp spin-density wave 
(SDW) order parameters which read as follows in terms of the Dirac fermions: 

P2kF = L'lRa, <S2kp = LiT^/sRp- (35) 

Using the representation (jCOp . we obtain the leading asymptotics of the 2kp CDW correlation 
function: 

(36) 
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where v and vq are respectively the spin- velocity in the Sp(2A'^) and Z^r sectors. We have obtained 
a similar estimate for the correlations which involve the 2kF SDW operator. At this point, these 
two correlation functions have the same power-law decay. The logarithmic corrections will lift this 
degeneracy and we expect that the SDW operator will be the dominant instability as in the = 1 



case 



381, 



39|. 



2. Confined phase 

Let us now consider the first spin-gapped phase, i.e. phase II of Fig. 1, which occurs when 
U <0 and V > NU/2 in the weak-coupling limit. In contrast to phase I, it has only one gapless 
mode which stems from the criticality of the charge degrees of freedom. The existence of a spin 
gap leads us to expect the emergence of a quasi-long-range BCS phase with the pairing of fermions, 



i.e a Luther-Emery liquid phase, as in the F = 1/2 case [38l . l39l |. However, this is not the case 
for F > 3/2 due to the existence of the Z^r symmetry ([5]) which remains unbroken in phase II. 
Indeed, it costs a finite energy gap to excite states that break this symmetry and the dominant 
instabilities must thus be neutral under Zjy. In particular, there is no dominant BCS instability 
in phase II since, as already stated in the introduction, the lattice singlet-pairing operator Pqq ■ 
([2]) is not invariant under the Z^r symmetry ([5]). Another way to see this is to use the low-energy 
description of the BCS operator obtained in Appendix C (Eq. (jClOj) ): 

^00 ~ Lij^^R^ ~: exp (i^27r/Ne,) : aiTr , (37) 



where Tr (p^^'^ takes a non-zero expectation value in phase II. Since the Z^v symmetry of the 
underlying Ising model is not broken, the Z^r Ising spin variable cJi has zero expectation value and 
short range correlations. Consequently, the equal-time correlation function of the BCS instability 
(1371) has a short-range behavior: {P^^ (x) Pqo (0)) ~ e"^'/?. The BCS singlet pairing is completely 
suppressed in this phase. 

In contrast, the dominant instabilities in phase II are expected to be: 

P2k, = LlRa, Mo = 6--^-i?t ^ . . . Ri^Ll ... (38) 

which are respectively the 2kF CDW and the uniform component of the lattice SU(2A)-singlet 
superconducting instability made of 2N fermions, i.e. the MS instability, Mj = ■..Cg^^. Both 
order parameters (j38|) are neutral under the Zat symmetry (j31|) and we notice that Mq is also 
invariant under the Z^r symmetry ()32p . The long-distance behavior of the correlation functions of 
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these operators can be determined using the identifications (|C9|) and (ICTSllCTl) of Appendix C: 

P2kp ~ : exp (i^J2^T/N^^ : /xiTV (39) 



N/2 

Mo ~ : exp (iV2^Q^ : J2 «P ^N/2-p Tr ^^^P), (40) 

where we have assumed, for the sake of simplicity, that N is even for the representation of the MS 
instabihty. The fields in the Sp(2A^) x Zat sector, which occur in these expressions, have non-zero 
expectation values in phase II so that: 



P2kp ~ :exp VW^^^c : (41) 



Mo ~ : exp (iV2TrNecj : . (42) 
We then deduce that these orders have power-law decaying equal-time correlation functions: 

(plkA^) P2kA0)) - x-^^/'', (43) 
(Mq^ (x) Mo (0)) ~ x-^/^^. (44) 

We thus see that CDW and MS instabilities compete and the key point of the analysis is the one 
which dominates. At issue is the value of the Luttinger parameter Kc- For Kc < N, the leading 
instability is p2kp which gives rise to a CDW phase whereas for Kc > N a MS phase is stabilized 
with a lattice order parameter Mj. Such a large value of Kc is not guaranteed for fermionic 
models with only short-range interactions as model ([3]). The full non-perturbative behavior of the 
Luttinger parameter as a function of the interactions U, V and density n is beyond scope of the 
low-energy approach. Its value should be determined numerically by means of Quantum Monte 
Carlo (QMC) technique or the density matrix renormalization group (DMRG) calculations. Such 
analysis has been performed recently for the spin-3/2 case, i.e. N = 2, and reveals that the MS 
phase (a quartetting phase in that case) exists for a wide range of attractive contact interactions 
at sufficiently small density [s^. In the general case, a strong-coupling investigation at small 



30[] . From the perturbative estimate (jl6p , 



density shows that an upper bound for Kc is -fCcmax = A" 
we see that Kc > 1 for attractive interactions and increases with \U\ and \V\ so that we expect 
the emergence of a MS phase in the general A^ case for sufficiently strong attractive interactions 
and small density. Such a phase is displayed on Fig. 2 and the doted line marks the crossover line 
between the CDW and MS phases. 

The MS phase with the formation of bound-states made of 2A^ fermions is remarkable and is a 
consequance of the higher-spin degeneracy when F > 3/2. The onset of this phase stems from the 
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existence of the Z^r symmetry which remains unbroken in phase II. This discrete symmetry leads 
to the suppression of the charge 2e BCS superconducting instabihty and confines the electronic 
charge to multiple of 2Ne i.e. the leading superfluid instability is a composite object made of 2N 
fermions. In this respect, the MS phase might be viewed as a kind of nematic Luther-Emery liquid 
since SU(2) invariant spin-1/2 electronic systems with a spin gap and gapless charge degrees of 
freedom, i.e. Luther-Emery liquids, exhibit a pairing phase 



3. Deconfined phase 

The properties of phase III are obtained from those of phase II with help of the duality symmetry 
"D (j27p on the fermions or (<I>c ^ Qc,'^N ^ '^n) in terms of the charge and spin degrees of 
freedom. Phase III has only one gapless charge mode and the Z^r symmetry, in the spin sector, 
is now spontaneously broken. The effective two-dimensional Ising model belongs to its low- 
temperature phase where the Z^v remains unbroken. We find now the emergence of a quasi-long- 
range BCS pairing phase which is described by the order parameter (|37p . Since cJi and Trcp^^^ 
acquire now a non-zero expectation value in phase III, we have the following leading behavior at 
low-energy: 

Pjo ~: exp (i^/2^^/N 6^) :, (45) 

which leads to the power-law decaying equal-time correlation function: 

(pJo(x)Poo(0))~x-V^^^^ (46) 

The Ztv symmetry, being spontaneously broken in phase III, does not confine anymore the electronic 
charge and thus accounts for the emergence of a BCS superfluid phase. In the spin-3/2 case, 



this phase has been found numerically by means of QMC and DMRG calculations [58t |. The 
MS instability, being neutral under Z^v, still has the power-law decay ([Hj) in phase III and is 
always subdominant with respect to the BCS operator (j45p . In contrast, the latter singlet pairing 
instability is now competing with the operator p^2Nkp '^^lich is obtained from the MS instability 
(f38l) by the duality transformation (f27|) : 

{27V) _ ai...ojv/3i---/3jv 7" K K t'^ t'^ ( A7\ 

t^2NkF ~ <^ai7i---c/aiv7iv -"'7i----"'7iV-^/3r--^/3iv' V^'/ 

which corresponds to a molecular density- wave (MDW) phase, with wave vector 2Nkp, made of 
2N fermions. Its low-energy description can be derived from Eq. (j40p by means of the Gaussian 
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duality ($c ^ ©c) and the KW transformation: 

N/2 

p'i^l ~: exp -.J^cip {-if^'-^er^/^-p Tr ^(^p), (43) 

p=0 

where we have used the KW transformation of the thermal operators of the Zjv CFT: ej {—lyej. 
The Sp(2A^) x Z^r fields in Eq. (j48p have non-zero expectation values in phase III so that the 
leading asymptotics of the MDW correlation reads as follows: 

{&ZH-)p?Zm-x---^. (49) 

Using Eq. (06]), we conclude that the BCS instability is the dominant one of phase III when 
Kc > 1/N whereas the MDW phase emerges for Kc < 1/N. The crossover line between these two 
phases is denoted by dotted lines in Fig. 2. 



4- Quantum phase transition 

As it has been discussed in section II, the nature of the quantum phase transition between the 
two spin-gapped phases II and III can be determined through the duality symmetry D (I27p . On the 
self-dual line g± = 0, i.e. 2V = NU, the low-energy properties of the transition are captured by the 
field theory (j28p . Using the conformal embedding approach, we observe that there is a separation 
of the Sp(2A^) and Zjv degrees of freedom in Eq. (j28p . Though the Sp(2A^) sector remains gapfull 
when U < 0, the effective Z^r Ising model is at its self-dual critical point and governs the phase 
transition. The Z^r quantum criticality might be revealed explicitly by considering the following 
ratio 7^ AT at equal times: 

TZn (x) = ((Pjo (^) ^00 (0))) /{a4 (x) Mo (0)). (50) 

In phase II, this ratio admits an exponential decay since the singlet-pairing instability is short 
range. Using Eqs. (011 H6|) . we find that TZj\f (x) ~ cte in phase III. At the Ztv quantum critical 
point, TZn (x) has a power-law behavior. Indeed, at this point, the BCS instability ([37]) simplifies 
as follows: 

Pjo (x) ~: exp (i^/2^^/Nec) : ai, (51) 

where we have averaged the Sp(2A^) degrees of freedom. Similarly, using the identification (|40p . 
we obtain that the MS instability Mq (x) is still given by Eq. (|42p . Therefore, we find that, at the 
critical point, TZn (x) has a power-law decay with a universal exponent: 

7^;v(x)~x-2^(^-l)/(^+2)^ (52) 
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i.e. an exponent 1 and 12/5 respectively for the Ising (N = 2) and three-state Potts {N = 3) cases. 
The long-distance behavior of the function 7^2 (x) close to the transition has been determined 
recently numerically by means of large scale DMRG calculations and the emergence of the Z2 
quantum criticality has been revealed 



58l |. However, the phase transition may be non-universal 
for larger N since we have neglected irrelevant perturbations in the weak-coupling limit |C/, <^ t 
that could become relevant at the quantum phase transition. Some insights might be gained 
from the symmetries of the problem. The relevant operators should belong to the Z^v sector, be 
neutral under the Z^v x Zjv symmetry and invariant under the KW transformation. The natural 
candidates are the thermal operators ej {j = 1, . . . ,N/2) of the Zjv CFT with scaling dimension 



+ l)/(-^ + 2) which transforms as ej {—lyej under the KW duality symmetry [361 ]. We 
thus deduce that, in the = 2, 3 cases, the quantum phase transition between phases II and III 
is universal and belongs to the Ising and three-state Potts universality classes respectively. For 
> 4, a strongly relevant perturbation £2 with scaling dimension 12/(A^ -|- 2) is generated in the 
Ztv sector. The resulting field theory which captures the quantum phase transition for F > 7/2 
becomes: 

Ses = Szff + X j (fx €2 (x) , (53) 

where stands for the action of the Z^r CFT. Model (I53p turns out to be an integrable de- 
formation of the Zat CFT 59|. The nature of the phase transition depends on the sign of the 
coupling constant A 5^. For A < 0, the field theories <\53h are massive and the phase transition 
is of first-order type. For A > it is known that model ()53p has a massless RG flow onto a 
Kosterlitz-Thouless (KT) U(l) gapless phase with central charge c = 1. 

Finally, we end this subsection by briefly considering the other phase transitions that occur in 
Fig. 2. The quantum phase transition between phases I and III is of KT type and its exact position 
isV = and f/ > 0, as it will be seen in section III B 4, which corresponds to the repulsive SU(2A^) 
Hubbard model. The transition between phases I and II is also expected to be of the KT type 
but its actual position is not clear. Within the weak-coupling approach of section II, its location 
is [/ = and ^ > but higher-order corrections might bend this line. 



B. Commensurate filling: one-atom per site 



We now discuss the nature of the phase diagram of model ([3]) for a commensurate 1/2N filling 
{kp = TT / {2N ao)) i.e. one atom per site. In this case, there is still a spin-charge separation ()13p 
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and the main modification occurs in the charge sector with the existence of an umklapp process 
which appears at higher order of the perturbation theory, and that transfers 2N fermions from one 
Fermi point to the other, i.e. l\RiL\R2 ■ ■ ■ L]2j^R2n + H.c. The low-energy Hamiltonian, which 
describes the charge degrees of freedom, becomes now a sine-Gordon model: 



Vr 



gc cos (VSTTiV ^>c) . (54) 



Such umklapp operator can also be determined from a symmetry analysis. Indeed, for the com- 
mensurate filling kp = 7r/(2A^ao), using the representation ([U]), the lattice one-step translation 
symmetry Tag reads as follows on the Dirac fermions: 

L„ ^ e-^"/2A^L«, Ra ^ e'""/^^ Ra, (55) 

so that we deduce ^ + \/vr/2A^ under 7^^ from Eq. ([9]) . The cosine term of Eq. (|54p is thus 
the one with the smallest scaling dimension and compatible with the translation symmetry. Since 
its scaling dimension is = 2NKc, we deduce that a charge gap is opened when Kc < l/N. A 
Mott transition occurs between a Luttinger phase and an insulating phase (see Fig. 3). In the 
latter phase, the charge bosonic field is pinned in one of the minima of the sine-Gordon model 



(^c) = \/^"^' 9c>0 



(^c> = ("^ + ^)' 9c <0, (56) 

m being integer. Three different Mott-insulating phases can then be defined depending on the spin 
degrees of freedom and the status of the Z^v discrete symmetry (see Fig. 3) . 

1. Mott I phase 

In the first Mott phase, the spin degrees of freedom remain gapless. This phase is qualitatively 
similar to the insulating phase of the repulsive SU(2A^) Hubbard chain at 1/(2A^) filling with 



2N — 1 gapless bosonic spin modes 



491 ]. In particular, it includes (when taking the limit of 



large repulsive U) the SU(2N) generalization of the Heisenberg model: 

n = J^Sfsf+„ (57) 

i 

where are the SU(2A^) spin operators on site i: = <^ai^ai3^P,i- This model is integrable by 
means of the Bethe ansatz approach and its critical properties are captured by the SU(2A^)i 
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FIG. 3: Phase diagram of model ([3]) for a filling of one atom per site. 



WZNW model 



48l |. The leading asymptotics of the spin-spin correlation functions in this phase 



can be easily obtained from Eq. (p6]) by gapping the charge degrees of freedom: 

-<Ar \ oB /n ^^\ COS {Ikpx) 



(5^ (x,r) 5« (0,0)) 



(58) 



2. Mott II phase 



The remaining Mott insulating phases of Fig. 3 are fully gapped. In region II, we have seen 
in section III that the Zjv (respectively Zjy) symmetry in the spin sector is unbroken (respectively 
spontaneously broken). We expect that the lattice translation symmetry Tag is spontaneously 
broken in the Mott II phase. This can be explicitly shown by considering the lattice 2A^-merization 
operator which is defined by: 



A 



2kp 



-iirj /N 



(59) 



In the continuum limit, this order parameter expresses directly in terms of the SU(2A^)i primary 
field: A2kp ~ T^g ~ iJ,iTT(j)^^\ Since the Z^r symmetry is broken, i.e. (/xi) / 0, Trg acquires a 
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non-zero expectation value in this phase: {/S.2kp) 7^ 0. We deduce, from the definition (j59p . that 
the translation symmetry is spontaneously broken and the Mott II phase displays a 2A^-merization 
with a 2N ground-state degeneracy. On top of this order, we also expect the formation of a long- 
range ordering of a 2kir CDW with order parameter p = cos{7rj/N) c]^jCaj- Using Eqs. (j6l4ip . 
we get the following identification in phase II: 

p ~ 4i?„ + RIl^ ~ : cos (v/27r/iV (60) 

which has a non-zero expectation value due to the pinning of the charge bosonic field (j56p . The 
2kF CDW coexists with the 2A^-merization operator (I59p in the Mott II phase. 



3. Mott III phase 

The Mott transition line Kc = 1/N coincides with the cross-over line between BCS and MDW 
phases in Fig. 2. For a filling of one atom per site, the BCS phase is still present while the 
MDW phase becomes a fully gapped Mott insulating phase. In this Mott phase, the Z^v symmetry 
is now spontaneously broken while the Z^v one remains unbroken. In the continuum limit, an 
order parameter of this phase is provided by the 2NkF MDW operator (j47p or the INkp CDW: 
P2Nkp = (^"'^'"'^'^^^'"^^ Laj^ . . . Llij^Rj3-^ . . . Rj3j^. In particular, using the character decomposition 
()C12p of Appendix C, we get: 

N/2 

P2NkF ~: exp (iV2^^c) ■■ «P ^N/2^p IV . (61) 

After averaging out the Sp(2iV)i x Zn fields, 2NkF MDW (08]) and 2NkF CDW order parameters 
have the same leading behavior which is given by: 



p2Nkp ~: exp (iV2^ $c : ■ (62) 



Using the position of the pinning of the charge bosonic field ([56l) . we deduce that (Re p2Nkp) / 
(respectively (Im p2Nkp) 0) when Qc > (respectively gc < 0). The electronic properties of this 
phase thus depends on the sign of the coupling constant Qc of the umklapp operator in Eq. 



This sign might be fixed using higher-order perturbation theory as in Ref. [Gl] or numerically by 
looking at the different lattice order parameters of the problem. In this respect, for gc > 0, we have 
the long-range ordering of a g = 2Nkp = vr/ao CDW: {p-,^) = (X^j q ("!)'' Q,Cj^Q,) 7^ while for 
< a spin-Peierls (bond) ordering is formed: (Osp) = {Ylj a(~l)''4+i a'^i>° ~'~ '^-'^■) 7^ 0. In any 
case, the Mott III phase is two-fold degenerate for all and the translation symmetry is broken. 
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This Mott phase is quite unusual since there is no one-particle charge density long-range fluctuation 
due to the spontaneously breaking of the Zjv symmetry. In the N = 2 case, at quarter-filling, it 
has been found numerically by means of the DMRG approach that a Mott transition occurs with 
the formation of a Akp = vr/oo bond ordering insulating phase [53]. We thus expect, at least for 
N = 2, that the coupling constant of the umklapp process is negative in phase III: gc < 0. Finally, 
we observe that the Mott III phase is the only fully gapped phase directly connected to the BCS 
phase (see Fig. 3). 

4-. Strong coupling expansion 

At large positive U, charge degrees of freedom are completely frozen by a large Mott gap, and 
it is possible to perform a strong coupling expansion that leads to an effective model for the low- 
energy spin degrees of freedom. Indeed, when t/U = in Eq. ([3|), the ground states with energy 
(the chemical potential is set to /i = U/2), which consist of states with all sites filled by exactly 
one particle ('singly-occupied states'), are well separated in energy from all other states by a large 
gap U. One should note that the gap to those states that are singly occupied on all sites except on 
two sites, one being empty and the other one being doubly occupied by a SU(2) singlet Pqq |0), is 
U + V. Therefore the following picture holds only when t <^ U,U + V to avoid formation of BCS 
pairs that occurs at large negative V. Nevertheless, we shall see that the BCS phase is captured 
by this approach, since it extends up to y = 0~. 

The huge degeneracy of the ground-state manifold consisting of singly occupied states is lifted by 
the hopping term 7io = —t ^ c^^ ^Ca,i + H.c. as soon as t/U departs from 0. To second order (the 
first non-trivial order) in perturbation theory, the effective Hamiltonian acting on singly occupied 
states reads: 

n^s = rno{i-r)-^ii-v)noP, (63) 

nu,v 

where V is the projector onto singly occupied states and Tiuy = 'H — Tio is the quartic part 
of the Hamiltonian The effective Hamiltonian can then be expressed in terms of SU(2A'^) 
spin operators defined on each site, Sf = c^ai^ap^P^i^ °^ which we can single out the Sf 
(a = 1, . . . , N{2N + 1)) that constitute a basis of the Sp(2A^) spin operators. One finds: 

4Ar2-l N{2N+1) 

^6/7 = ^0^ ^ SfS^^i + Gi'^ ^ SiSi_^i, (64) 

i A=l i a=l 

with the couphngs Go = f + cttf) and d = |J - ^) . 
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When V = 0, Gi = and the model reduces to the SU(2A^) antiferromagnet in the fundamental 
representation. As recalled in section IIIIB 11 this integrable model is described at low energy by 
the SU(2A^)i WZNW model perturbed by a marginally irrelevant current-current interaction. At 
small V/U, i.e. in the vicinity of the SU(2A^) invariant point, we can thus obtain a continuum 
limit of Hamiltonian ()64p in the form of a perturbed SU(2A^)i WZNW model. The calculation 
relies on the continuous representation of the spin operators, 5^ ~ S^{x = jao) = II ~^ -^R ~^ 
Yl'^^i^ '^nie^*™''^^^5^^^, where the numbers Om are non-universal amplitudes. 

This procedure yields a low-energy Hamiltonian density exactly of the form (|20p , with couplings 
^11 ([/, V) and g±{U, V). The actual value of g^^ and g± are not known (they depend on non-universal 
amplitudes), but it is sufficient to know that g^^(U,0) = g±{U,0) < at the SU(2A^) invariant 
point, and that (jy — oc Gi oc ^ at small V/U. The RG analysis of section Hi Bl thus leads to the 
conclusion that the spin sector undergoes a transition at V = 0, between a quantum critical phase 
(y > 0), that coincides with the Mott phase I, and a spin gapped phase {V < 0) with unbroken 
symmetry, which is Mott phase III. The strong coupling expansion thus establishes that the 
exact position of the transition line between these two Mott phases is U > 0,V = 0. 

IV. CONCLUDING REMARKS 

In summary, we have studied the zero-temperature phase diagram of one-dimensional spin- 
F = N — 1/2 fermionic cold atoms with contact interactions. In the low-energy limit, a CFT 
approach has been developed to deduce the main physical properties of the model by exploiting 
the presence of an extended Sp(2A^) symmetry. In the F = 3/2 case, this Sp(4) ~ S0(5) symmetry 
is exact on the lattice while for larger F its existence stems from a fine-tuning of the coupling 
constants i.e. the scattering lengths of the atoms. 

The phase diagram of the model is very rich for incommensurate filling with the competition 
between instabilities of very different nature. In particular, two different superfluid phases are 
found for attractive interactions: a Sp(2A^)-singlet BCS pairing phase and a MS phase. The latter 
phase is an SU(2A^) singlet which is formed from bound states of 2A^ fermions. For instance, for 
N = 2, i.e. F = 3/2, it is the analogous of the a particle in nuclear physics. At the heart of 
this competition is a Z^r discrete symmetry which is the coset between the center of the SU(2A^) 
group and the center of the Sp(2A^) one. Since Sp(2A^) and SU(2) share the same center group, 
this Ztv symmetry is not an artifact of the extended Sp(2A^) symmetry for F > 7/2 and is present 
for general SU(2)-invariant multicomponent Fermi gas. The Zjy symmetry plays a crucial role and 
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provides a new important ingredient which is not present for F = 1/2. 

In the low-energy approach, this discrete symmetry corresponds to the symmetry of the low- 
temperature phase of an effective two-dimensional Zjv generalized Ising model. If the Z^r symmetry 
is not spontaneously broken, we find the suppression of the singlet-pairing BCS phase with the 
confinement of the electronic charge e to multiple of 2Ne. In this case, the nature of the dominant 
instability of this phase depends on the non-universal Luttinger parameter K^. If < N, the 
leading instability is a 2kp CDW while a quasi-long-range MS phase emerges for Kc > N and 
should be the generic phase for attractive interactions at sufficiently small density. Due to its 
confinement properties and the existence of a spin gap, this MS phase can also be viewed as nematic 
Luther-Emery liquid. When the Z^r symmetry is spontaneously broken, a quasi- long-range BCS 
phase is stabilized when Kc > 1/N whereas a MDW phase appears for Kc < 1/N. The quantum 
phase transition between the confined and BCS phases is described by the Z^r parafermionic CFT 
perturbed by the second thermal operator. For F = 3/2 and 5/2, the transition is universal 
and belongs respectively to the Ising and three-state Potts universality classes. For higher spins 
F > 7/2, the transition is non- universal and is either of first-order or in the U(l) KT universality 
class. 

For commensurate filling of one atom per site, in sharp contrast to the F = 1/2 case, a Mott 
transition is predicted to occur for > 3/2 as the result of the spin degeneracy. Three different 
Mott-insulating phases are expected which depend on the status of the Z^v symmetry. A Mott 
phase with 2N — 1 gapless spin modes is first found which is qualitatively similar to the insulating 
phase of the SU(2A^) Hubbard chain. The two others have a spin gap and display spin-Peierls 
(bond) ordering with a 2N or a two-fold ground-state degeneracy. 

Regarding perspectives, it will be interesting to investigate numerically the phase diagram for 



F = 5/2 spins by means of large-scale DMRG or QMC calculations [62l |. In particular, the three- 
state Potts universal behavior of the quantum phase transition might be confirmed numerically 
as well as the study of the Mott transition and the different Mott-insulating phases for one atom 
per site. From the theoretical point of view, another important question is the phase diagram for 
incommensurate filling of the SU(2) model ([T]) for > 3 without any Sp(2A^) fine-tuning. Finally, 
we hope that the exotic MS phase for general F = N — 1/2 spins, discussed in this paper, will be 
observed in future experiments in ultracold spinor fermionic atoms. 
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APPENDIX A: GROUP THEORY CONVENTIONS 

The aim of this Appendix is to define the algebraic conventions used in the low-energy approach 
as well as the RG calculations which are presented in Appendix B. 

Let us start with the SU(2A^) group. The generators of the Lie algebra of SU(2A^) in the 
fundamental representation (2A^ x 2A^ matrices) are denoted by = 1; • • • j 4A^^ — 1) with the 

normalization: Tr(r^T^) = 6^^ /2. They satisfy the commutation relation: 

[r^,r^] =i/^^^r^, (ai) 

jABC ]-)gi]2g the structure constants of SU(2A^). An useful group identity for the derivation of the 
continuum limit of the lattice model ^ is: 

The SU(2A^) generators can be divided in generators of the subalgebra Sp(2A^), (a = 
1, . . . , N{2N + 1)), and generators of the complement of Sp(2iV) in SU(2iV): (i = 1, . . . , 2N'^ - 
N - 1). They are normalized as: Tr(r"r'') = S""^/!, Tt{T^T^) = 6'^ /2. The various structure 
constants are then given by: 

which defines a Z2 graduation of the SU(2A^) Lie algebra. These coefficients can be determined by 
expressing the generators T"^, in terms of a direct product between SU(A^) and SU(2) generators. 
The SU(A^) generators admit a simple NxN matrix representation which falls into three categories: 

• Symmetric part: 

=l{6,aS,f3 + SipSja) {l<i<j<N) (A4) 

• Antisymmetric part: 

{mS^) = -^('5-'^./3 - S^(36,a) il<i<j<N) (A5) 



= if^^T^, [T\ T^] = if^^T", [T", r] = if'^T^, (A3) 



28 



• Cartan generators: 

i^^)a0 = } I Tl ^ak^Pk - m 5a,m+l^p,m+l ) , (m = 1, . . . , iV - 1). (A6) 

V2m(m + 1) ] 

The Sp(2A^) generators T"" (a = 1, . . . , N{2N + 1)) are then given by the following set: 

= I ^ iv4^^ <?, ^ M^f 0/2, ^ a, ^^/iv a I , (A7) 



being the Pauli matrices. The metric Jap^ which defines the Sp(2A^) group, takes a simple form 
in this scheme: J = In® {—i<^2)-, and is thus the generalization of the antisymmetric tensor eap for 
higher spin F > 1/2. Using the representation (|A7p . we also obtain the following identity which is 
useful for the derivation of the low-energy approach of section II: 

^a/3^7P = ^ i^apSfB^y - Ja^Jf3p) ■ (A8) 

a 

Finally, the remaining generators T* (i = 1, . . . , 2A^^ — — 1) read as follows: 

T^ = \^^M^f^a, -Lm«®I2, ±M^^hY (A9) 

APPENDIX B: TWO-LOOP RG RESULTS 

In this Appendix, we present the two-loop RG equations of the current-current model (I20p : 



n 



AttJ 



^j9±9\\ Ng± / 2 ^ 2 , 
9± = - ( <7|j + 5± ) , 



(Bl) 



2tt (47r) 

where y = dg±^^^/dt, t being the RG time and we have neglected the spin-velocity anisotropy. To 
establish this result, we use the notations of Ref. (we modify normalizations for our purpose) 
and define the operators: 



in terms of which the interacting part of the action reads 



Si = (fx QaOa- 



(B2) 



(B3) 
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Structure constants C, D, C are defined for the algebra of the local operators Oq, and associated 
pseudo stress-energy tensors =: I^^^I^l '■ (-2) and T±{z) =: '■ (z), according to the 

OPE's: 

O^{z,z)Of,{0,0) ~ -l--C2.O^(0,0), (B4) 



T„(z)O^(0,0) ~ -i-^ (2kDl, + Cl,)o^{0,0). (B5) 



87r2 |z|2 ^"/3^7\ 

Here, to disentangle the structure constants D and C, we need to vary fictitiously the level k of the 
Kac-Moody algebras (A; = 1 in our case), thus considering the more general problem of perturbing 
a SU(2A^)fc model by marginal current-current interactions built on the currents of the Sp(2A^)fc 
that are embedded in it. We still write for the S\J{2N)k currents. The only formal difference 
appears in the current-current OPE's: 

u ^AB A fABC 

with a similar expression for right currents. 

The two-loop /9-function takes a compact form in terms of the structure constants 

5a = ^ gpg, - ^ C^c^D^; 9,9,9u. (B7) 

The structure constants are calculated to be (we give only the non vanishing entries): 

Cj'll = -2(iV + l), C7||'^ = -2(iV- 1), C7]_'^ = C^'ll = -2iV 
Cj'" = 2(A^ + 1), C\_'^ = 2N, cl^ = 2{N-l), cf'^ = 2iV, (B8) 



so that Eq. (|B7p gives the two-loop result p3T 



APPENDIX C: Zat PARAFERMIONIC CFT 

In this Appendix, we discuss the Z^v parafermionic CFT and present the technical details of the 
conformal embedding approach of section II which are used in section III for the determination of 
the phase diagram of model 

The parafermionic CFT describes the critical properties of two-dimensional Zat generalization 
of the Ising model. The lattice spin ar of such a model takes values: e*^™/''^, m = 0, . . . , A^ — 1 and 
the corresponding Hamiltonian is Zjv invariant. Using the KW duality symmetry, the spins are 
replaced by the dual spins /Uf which belong to the dual lattice and they define a Z^v symmetry. Two- 
dimensional Zat Ising models display two gapped phases, i.e. the ordered (respectively disordered) 
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phase where the Ztv (respectively Zjv) symmetry is spontaneously broken, as well as a quasi- 
long-range KT phase for prime N > 5. The critical properties, along the self-dual manifold, are 
captured by the Z^v x Zjv parafermionic CFT with central charge c = 2(A^ — 1)/(A^ -|- 2) [36,]. In 
the scaling limit, the conformal fields and fik = cTAr_jfc, = ^N-k) with scaling dimensions 
dk = k{N — k)/N{N + 2) {k = 0, . . . , N — 1) describe the long-distance correlations of a^. and jj^. 
In particular, they carry respectively a (fc, 0) and (0, k) charge under the Z^v x Zjv symmetry: 

^i2nmk/N^^ Under Zn 
^ e^^-^'^/^^fc under Zn, (CI) 

with m = 0, . . . ,iV — 1 and ak (respectively fik) remains unchanged under the Zjv (respectively 
Ztv) symmetry. 

The Ztv parafermionic CFT is generated by the parafermionic currents "^kL a-nd "ifkR (^II b. ~ 
^N~kL,R) which are holomorphic and antiholomorphic fields with conformal weights: = hk = 
k{N — k)/N. They carry respectively a {k,k) and {k,—k) charge under the Zjv x Z^v symmetry 
and they transform as follows under the KW duality symmetry: 

"fkL ^ ^-fcL, ffci? ^ (-1)' ^L, (C2) 

n addition, these currents satisfy the parafermionic algebra which is defined by the following OPEs 



"fkL (z) -^gL (w) ~ Ck,g {z - c^)-2^'?/^ ^k+qL (c^) , for A: + g < 4 (C3) 

^kL (Z) (^) ~ Ck,N-q {Z - uy^li^-^y^ ^k-qL (^) (C4) 

^kL (z) (-) ~ - (l + ^ - ^ (-)) ' (C5) 

with similar results for the right parafermionic currents. In Eq. (jCSp . T(z) denotes the stress- 
energy tensor of the Zjv CFT and the numerical coefficients are given by: 

2 r(fc + g + l)r(5-fc)r(5-g) 

T{k + l)T{q+l)T{6-k-q)T{5)' ^ ' 

Apart from these fields, the Zjv CFT also contains neutral fields, i.e. invariant under the Z^v x Z^v 
symmetry, ej {j = 1, . . . , [N/2]) with scaling dimension Dj = 2j{j + 1)/(A^ + 2) which are the 
thermal operators of the theory. 

It is well known that the Z^v CFT can be viewed as a coset (see for instance Refs. [s^,!^^): Ztv ~ 
SU(2)Af/U(l)Ar or Zn ~ [SU(A^)i x SU(A^)i]/SU(A^)2. In fact, there is another coset description 
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of the Zjv CFT which is crucial for the low-energy approach of section II: Zjv ~ SU(2A^)i/Sp(2A^)i. 
The corresponding character decomposition has been found by Altschuler [56| and reads as follows: 

N N 

X (A,) = B^X (A,) =Yv4x (A,) , (C7) 

3=0 3=0 

where x i-^k) = 0, . . . , 2N — 1 (respectively x (Aj)) are the characters of the conformal tower 
corresponding to the SU(2A^)i (respectively Sp(2A^)i) primary fields ^^''^ (respectively (f)^^^) with 
scaling dimensions: k{2N-k)/2N (respectively j{2N + 2- j)/2{N + 2)). In Eq. jCTl), are the 
branching functions of the coset which identify to the level- string functions of the SU(2) at 



current algebra 



6j]. The latter are related to the branching functions fr^ of the coset Z 



N 



SU(2)jv/U(l)iv: b^f. = {r] being the Dedekind function) with the constraint j = k (2) 64|. 

The identification ()C7|1 and the non-Abelian bosonization rule ([9]) enable one, in principle, to 
express any fermionic operators of the U(2A)i CFT in terms of a U(l) charge boson field and 
operators of the Sp(2A)i x Zjv CFT. Let us first consider the example of the representation of the 
U(2A)i field: llRa. Using Eq. ([9]), we have: 



LlRa ~: exp (i^27r/N^A : TV g. (C8) 



The next step of the approach is to use the decomposition (\C7\i with k = I. The Sp(2A)i primary 
field (f>^^\ which transforms in the fundamental representation (self-conjugate) of Sp(2A), will 
appear as well as a Z^r field with scaling dimension (N — 1)/N{N + 2) i.e. ai or /ii. One way to 
resolve this ambiguity is to use the parity symmetry P in the continuum limit. Under the parity, 
{L,R)a — > {R,L)a and ^cL,r — ^ —^cR,l so that Tt g ^ Tr g"^ from Eq. (|C8p . However, it is 
known that, under P, Uk remains invariant while fik fJ-^ Isfil- We thus deduce the following 
correspondence: 

Ll^Ra ~: exp (^iv^27r/A$c) : Tr g ~: exp (iy^2^<^c) ■ /^iTr (p^^\ (C9) 

Another interesting operator is the Sp(2A^) singlet LaJ'ajjR^^ field whose representation can be 
obtained from Eq. (jCOp with help of the duality symmetry ()27l) of section II. Since the latter 
symmetry corresponds to the Gaussian duality and the KW transformation of the Z^r Ising model, 
i.e. (Tfc <-> /Uj;, we have: 



LlJafsRl ~: exp [iy^2^ec) : aiTr cj)^^^ , (CIO) 
which is consistent with the parity symmetry since LliJapR^p is invariant under P. 
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A more complicate example is the identification of the operator Mq = 
e^i-^'v^i-^^i^ii ...fii^Lj;^ ...lJ;^ in the U(l) x Sp(2iV)i x basis. This field is an 
SU(2A^) singlet and plays a crucial role in the discussion of the physical properties of the spin-F 
fermionic model ([3|) since it is the continuum representation of the MS instability. Using the 
non-Abelian bosonization rules Q, we first find that: 

Mo ~: exp (^iV27rA^ec) : Tr <^^^\ (Cll) 

^^^^ being the SU(2A^)i primary field which transforms into the self-conjugate representation of 
SU(2A^) given by a Young tableau with N boxes and one column. The next step of the approach 
is to use Eq. ()C7p with k = N. Let us assume that is even, i.e. A^ = 2n, so that we find in this 



6J. It 



case: 

n 

X(A2„)=??5^C^Px(A2(n-p)), (C12) 
p=0 

where we have used the symmetry property of the level-A" string functions: = c^Z^ 
is known that rjc^^ with p = 0, . . . , A^/2 describe the conformal towers corresponding to the 
thermal operators ep of the Z^v parafermionic CFT |6^] . We thus deduce the following identification 
of the MS instability for even A^: 

N/2 

Mo ~: exp (iV2^A^ec) : ^ ap ejv/2-p Tr ^(^p) ^ (C13) 

Up being some numerical coefficients which are not important in what follows. When A^ is odd, a 
similar approach gives: 

iN-l)/2 

Mo ~: exp (^^2^6,) : ^ e(Ar_i)/2-p Tr (/^(^P+i). (CM) 

p=0 

Finally, we end this Appendix by mentioning an interesting identification of the first 
parafermionic current ^'i_r,l which generates the algebra (jC3l IC41 IC5|1 . First, we observe that 
LaJ^afsL^ls and R^iJapR^ are chiral Sp(2A^) singlets so that only the identity operator in the Sp(2A^) 
sector can occur in their representation. We then find the following identifications: 

vr 



LiJapL\ ~ ^ : exp (i^Sn/N ^^l) : ^IL 



RijapR\ ^ ^ • {-W^T^/N ^cr) ■■ ^IR, (C15) 
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which satisfy the parafermionic algebra (jCSP for k = 1. The representation of the parafermionic 
current ^kR,L with k > 2 can then also been derived from the correspondence (|C15h and the 
application of the parafermionic algebra ()C3|C4p . 
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